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SUMMARY 


A real-time simulation of the UH-60A Black Hawk Helicopter based on the 
Sikorsky Gen Hel mathematical model is compared with flight test data and Sikorsky's 
nonreal-time computer program. An overview of the mathematical model is given, 
including a description of updates and expansion of the model by NASA. Also, real- 
time programming requirements and techniques are explained. The simulation displays 
very good agreement with flight test data and excellent agreement with the nonreal- 
time program on which it is based. The large dedication of effort required to 
develop the real-time program suggests a need for development of generalized tech- 
niques for real-time simulation of high-fidelity physical models of rotorcraft. 


INTRODUCTION 


High fidelity modeling of helicopter dynamics is increasingly becoming a 
requirement for pilot-in-the-loop handling qualities research. Dynamics once con- 
sidered to be of secondary importance, such as vehicle response to rotor speed 
variations, are now seen to be areas where improvement is desirable for future 
generations of rotorcraft. In addition, high-gain feedback control systems require 
considerable model sophistication. If careful consideration is not given to high- 
order dynamics, possible instabilities of the actual closed-loop system will not be 
reproduced. Simulation models are therefore required which provide the fidelity 
necessary to reproduce higher-order effects over a large flight envelope in real 
time. 


High-fidelity engineering simulations which are based on physical representa- 
tions of dynamic systems have advantages over linear-map or parametric equation 
representations for many simulation investigations. In this type of model, internal 
physical quantities such as individual rotor blade positions or shear forces at a 
hinge are determined. The dynamics of each helicopter subsystem and its effects on 
other subsystems are modeled in physical terms. This allows a user to modify a 
component subsystem, such as a lag damper, based on tests of the actual hardware. 

If an interaction between two or more subsystems has an effect on handling qualities 
(for example, the torque-balance interaction of the main rotor, tail rotor, engines, 
and fuel control system), the high-fidelity physical model may be used to investi- 
gate this effect over a wide range of flight conditions. Simpler models may cor- 
rectly simulate the interaction only about a narrow range of conditions for which 
estimates of each interacting quantity are valid. High-fidelity piloted simulation 
may also be used to investigate handling qualities or controllability of the modeled 
aircraft outside the normal range of operation. Because vehicle systems are modeled 
physically, internal quantities can be observed which may be important for reasons 
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not related to handling qualities, such as flight safety. For example, the proxim- 
ity of each blade tip to the tail boom may be monitored as the pilot operates in an 
off-design condition. 

A high-fidelity real-time model of the Sikorsky UH-60A Black Hawk helicopter 
has been under development for several years at NASA Ames Research Center for use on 
the Vertical Motion Simulator (VMS). Validation of the mathematical model of the 
helicopter was initially performed by Sikorsky Aircraft (ref. 1). Validation of the 
real-time model is necessary to establish that real-time computer coding considera- 
tions do not dynamically alter the model and that subsequent updates to the model 
have in fact further improved the fidelity of the model. 

The real-time simulation is first demonstrated to-be statically valid over a 
typical envelope of operation. Examples of dynamic response to pilot input are then 
used to show concurrence with the previous validation study and to illustrate the 
effects of improvements by comparison with flight-test data. The work presented 
here should not be interpreted as a complete validation over the entire operating 
range of the simulation. 


MATHEMATICAL MODEL DESCRIPTION 


The real-time mathematical model of the UH-60A is based on the generalized, 
modularized programs developed by the manufacturer for time-domain analysis of total 
helicopter systems. These programs make up the Sikorsky General Helicopter Flight 
Dynamics Simulation (Gen Hel). A real-time version of this mathematical model was 
provided to NASA Ames under contract (ref. 2). From this model, a computer program 
was written to execute in real time on the CDC 7600 computer. In addition, the 
fidelity of the model was increased as needed, notably in the area of engine and 
drive train modeling. 

The model is a nonlinear, representation of a single main rotor helicopter that 
is accurate for a full range of angles of attack, sideslip, and rotor inflow. Six 
rigid-body degrees of freedom as well as the main rotor flapping, lagging, air-mass, 
and hub-rotational-speed degrees of freedom are modeled. Because of its modular 
nature, each major force- and moment-producing component of the aircraft is treated 
as an independent entity. The underlying framework of the program is the interfac- 
ing between these components. All of the interfaces are physical quantities such as 
forces, moments, attitudes, body-fixed velocities, and downwash velocities. In this 
way, interaction between the components is modeled while maintaining the flexibility 
to modify or replace any one component without modification to the other compo- 
nents. Figure 1 is a block diagram representation of the simulation components. 

A blade-element approach is used to model each main rotor blade. Total rotor 
forces and moments are produced by summations of forces from each blade, which are 
determined from aerodynamic, inertial, and gravitational forces. Aerodynamic forces 
are computed from angle of attack and dynamic pressure acting on each blade segment 
based on the orthogonal velocity components. These components are determined as 
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Figure 1.- Gen Hel simulation components. 


functions of blade azimuth, lag and flap angles, local velocity of the blade seg- 
ment, and local downwash. Downwash is approximated to have a first harmonic distri- 
bution as a function of wake skew angle. Blade inertial and gravitational forces 
are computed from blade rotational velocity, lagging and flapping velocities and 
accelerations, and blade position. Wo dynamic twisting or bending of the blades is 
modeled, although a preformed blade twist is represented through adjustment of the 
geometric pitch of each segment. The summations of forces act on the airframe at 
the blade hinge and lag damper locations. Rotor moments result from blade hinge and 
lag damper offsets from the main rotor shaft. 
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The tail rotor thrust is represented by linearized Bailey theory. Interference 
effects from the aerodynamic modules are accounted for as empirically determined 
blockage factors. The main rotor downwash is used to modify the tail rotor inflow. 
The aerodynamics of the fuselage, stabilator, and vertical tail pylon are each 
represented in separate modules so that the nonlinear interference effects of the 
main rotor and the interference between the components are modeled separately. The 
aerodynamic function tables were developed from wind-tunnel test data, and refer- 
ence 3 was used to extrapolate and modify the available data. 


MODEL UPDATES 


Some inertial terms of the blade flapping acceleration equation had been 
dropped by Sikorsky under simplifying assumptions. During real-time implementation 
of the RSRA mathematical model (ref. 4), it was found to be necessary to complete 
the inertial terms. The Gen Hel RSRA blade-element rotor model is identical in 
formulation to the UH-60A model; all additional terms were therefore incorporated in 
the UH-60A simulation. A complete description of the added flapping-equation iner- 
tial terms may be found in reference 5. 

As part of an ongoing VMS effort to explore the improvements attainable in 
helicopter handling qualities by improving rotor speed control, a real-time digital 
simulation for the GE T700 engines and their hydromechanical and electrical fuel 
control systems was developed by NASA Ames. Based on a simplified hybrid computer 
program developed by NASA Lewis, it was expanded to a level of sophistication neces- 
sary for compatibility with the Sikorsky Gen Hel mathematical model. Rather than 
using a partial derivative representation (as was done in ref. 2), the major engine 
components and dynamics of several internal engine states are modeled on the basis 
of physical laws relating to transfer of energy and mass flow. The simulation is a 
close representation of T700 performance and dynamics over the full operating range 
and its component formulation allows great flexibility. 

Because of this increased level of detail, it was necessary to incorporate a 
modified version of the expanded gearbox model developed by Sikorsky for the RSRA. 

It allows an accurate representation of power dissipation by the main rotor, tail 
rotor, accessories, and transmission. The inertias of all of the major rotating 
components are also modeled. The tail rotor module was expanded to include a 
torque-required computation by the addition of Bailey drag terms (described in 
ref. 6). A new formulation of the engine-overrunning clutch was also developed by 
NASA Ames to allow correct operation at flight idle and in situations involving loss 
of engine power. 

In early real-time qualitative assessments, the pilots commented about unreal- 
istic lags between collective control input and vehicle response and their inability 
to perform precision hover tasks without causing pilot-induced oscillation. After 
examination of total-system and model delays and lags, it was determined that the 
bandwidth of the transfer function representing primary servo dynamics was too 
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low. The swashplate angle of the actual aircraft is controlled by three primary 
servos which exhibit nonlinear characteristics, both individually and when working 
together. However, the linear transfer function representation is considered to be 
a level of sophistication which is appropriate for use with Gen Hel. After consul- 
tation with the manufacturer, a bandwidth which is based on the no-load acceptance 
testing of the actuator was used in place of the original math model implementa- 
tion. As will be shown, the time histories of the transient responses reflect the 
improvement. Pilot comments also suggest that collective response characteristics 
are much improved. 

Appendix A contains the complete equations describing the model updates imple- 
mented on the Ames real-time simulation. Variables used in the appendix are 
Sikorsky Gen Hel variable names (ref. 2). 


REAL-TIME CONSIDERATIONS 


Real-time simulation of a blade-element rotor requires that time steps be as 
small as possible. Modeling of high bandwidth dynamics also requires small time 
steps. In addition, delays from cockpit input to visual- and motion-system output 
must be minimized. The outputs from each program module must be computed from one 
pass through the module. Internal iteration is used only when it is unavoidable, 
because a constant time step must be based on the maximum number of iterations 
needed. For each program pass, forces and moments from each component are computed 
and summed, from which the net translational and angular accelerations acting on the 
airframe are determined. The resulting body-fixed velocities and positions are used 
sequentially in the succeeding intervals to emulate a continuum solution to the 
total system. 

The time-step size has a large effect on the output from the blade-element 
rotor program. The rotor model contains its own integration algorithms, tuned to 
give correct flapping and lagging positions and velocities in the rotating reference 
frame of the rotor hub. Because each rotor blade is modeled at its azimuth orienta- 
tion, each update results in an azimuthal transition which, if too great, causes 
integration errors within the rotor that result in contaminated rotor outputs. 

Also, as in actual rotors, the total rotor outputs contain high frequency harmonics 
with a fundamental frequency equal to the number of rotor blades times the angular 
velocity of the hub (referred to as the N/rev frequency). A maximum step size must 
therefore be set to prevent aliasing of the fundamental frequency into lower 
frequencies . 

The sequential nature of the "one-pass requirement" leads to algebraic loops in 
digital simulation. Algebraic loops are a result of sequential computation of 
interdependent parameters and must be prevented wherever possible through careful 
designing of code sequence and hierarchy, A discrete model must necessarily provide 
a value for a variable at the end of a time interval based on both its previous 
value and the values of parameters upon which it is dependent at the beginning of 
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the interval. In the real world, the variable is dependent on the values of other 
parameters at the same instant in time. The occurrence of an algebraic loop can 
result in significant differences in phase compared to the continuous system and the 
dynamic stability of the system can be affected. The real-time program has been 
sequenced, integration algorithms have been chosen, and the time-index corrections 
have been used to minimize these effects. 

For example, the rotor speed degree-of-freedom results in an algebraic loop 
between the rotor, engine, and fuel control system which causes a significant propa- 
gation of error with increasing time step. Changes of rotor speed are a result of 
an imbalance of the main rotor torque required and the engine torque available. The 
hub angular acceleration is integrated to obtain hub velocity for each pass in the 
gearbox module, which is sequentially executed after the rotor and engine modules. 

At the beginning of each pass through the rotor, blade lag velocities and positions 
are computed from the integration algorithms. Lag acceleration, which is function- 
ally dependent on hub acceleration, is then computed. However, the hub acceleration 
is actually an output from a previous pass through the gearbox and is therefore 
delayed by one interval. As explained in reference 7, rotor torque harmonic varia- 
tions are largely a function of blade lag angles. Because rotor torque is used in 
the gearbox to determine the value of hub acceleration, an algebraic loop is 
formed. Real-time coding therefore requires that lag accelerations be corrected 
before their use in the blade-lag integration algorithms. The latest value of hub 
acceleration, calculated after the previous rotor pass, is used to modify the lag 
accelerations, and the effects of the previous hub acceleration used are subtracted 
out. This results in sequentially correct values of lag angle and resultant 
torque. Figure 2 illustrates the effect of this change for a time step of 
20 msec. An abrupt decrease of collective pitch at 60 knots causes a variation in 
rotor hub speed and acceleration as shown in figure 2(a). Hub acceleration oscil- 
lates at the M/rev frequency, corresponding to the physical changes of torque at the 
hub. Without the correction, the algebraic loop causes a spurious 10-Hz oscillation 
at a large amplitude (fig. 2(b)). 

Implementation of the blade-element rotor model also led to a reexamination of 
the basic kinematic model used at NASA Ames for all piloted simulation (ref. 8). 

Body translational and angular accelerations are computed by this module, but they 
are not normally updated for synchronization with the following time interval 
whereas velocities and positions are. Because the rotor module must use body accel- 
erations as well as velocities and positions to determine hub accelerations, an 
algebraic loop is generated. It was found that sensitivity of the rotor module to 
errors in phase is greater than its sensitivity to small magnitude errors. The 
kinematic model was therefore modified to provide extrapolated accelerations, 
thereby eliminating much contamination of rotor output. 

The "mixed time indices" which are the inevitable result of sequential one-pass 
coding have been minimized for the entire program using methods as described 
above. Correction of mixed indices alleviates the need for filtering or special 
numerical techniques for attenuation of rotor noise as were used in McFarland 
(ref. 9). The program can be executed with time steps as great as 20 msec with no 
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Figure 2.- Rotor hub speed and acceleration response to a 2-inch down collective 
step at 60 knots (0.02 sec time step), (a) With blade lag acceleration computa- 
tion sequence correction; (b) without correction. 






measurable effect on blade dynamics. Temporally mixed index values that still exist 
were determined to be too computationally costly to correct, but their influence on 
system dynamics is minimal. The real-time program sequence and hierarchy are shown 
in figure 3. 



Figure 3.- Real-time program sequence and hierarchy. 
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FLIGHT-TEST DATA 


The test data used for comparison were obtained in a series of tests conducted 
by USAAEFA for use in validation of the Rotorcraft Systems Integration Simulator. A 
detailed description of the testing is given in reference 10. The trim data and 
selected transient-response time histories were provided to Sikorsky for use in 
their validation of the mathematical model (ref. 1). The test aircraft, 

S/N 77-22716, was an early production Black Hawk containing USAAEFA test instrumen- 
tation. An airspeed boom was mounted forward of the nose, and a movable ballast 
cart was used to adjust the longitudinal center of gravity as needed. In addition, 
blade pitch, flap and lag instrumentation was provided by Sikorsky. 

Data were obtained while maintaining a constant thrust coefficient at specified 
vehicle and blade-tip Mach numbers. The altitude was increased as fuel was expended 
and the rotor speed was varied as a function of temperature. Fuel crossfeeding and 
the ballast cart were used to maintain constant lateral and longitudinal c.g.'s, 
while the waterline c.g. was allowed to vary. 

Wo stability augmentation was used during transient-response data runs. Analog 
and digital stability augmentation systems, the flightpath stabilization system, and 
the stabilator control system were disabled and the pitch bias actuator was centered 
and locked. The stabilator was fixed for each run; its position was determined as 
the stabilator control system would have set it for the aim airspeed and collective 
determined by the specified thrust coefficient. This is a highly degraded configur- 
ation; the results are not representative of a UH-60A in normal operation. The 
configuration, however, allows a validation of the basic vehicle model without 
masking model deficiencies, as may occur when control system models are operating. 

The transient responses consisted of individual -axis steps, pulses, or doublets 
of one inch or less in both directions. The velocities varied from hover to 
140 knots true airspeed. The test procedure normally consisted of stabilizing in 
trim with one of the two redundant stability augmentation systems on; this was 
disabled one second before the control input. Unsatisfactory stability characteris- 
tics of the unaugmented aircraft, especially in pitch, required the pilot to ini- 
tiate recovery within a few seconds of the input. 


STATIC TRIM 


The steady-state trim characteristics of the two simulations and the flight 
test data are given in appendix B. A general six degree-of-freedom trimming algo- 
rithm was used to determine iteratively the steady-state flight conditions by 
adjusting aircraft attitude and control positions to zero all net forces and moments 
acting on the aircraft. The rigid-body equations of motion were bypassed during 
trim, but the rotor module was allowed to step forward in time. The aircraft con- 
figuration was set up to duplicate static-trim cases performed by Sikorsky to match 
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USAAEFA test data. Flight test values of aircraft weight, c.g. location, and trim 
main rotor speed were used. Comparisons of flight test values of stabilator inci- 
dence angle with known values based on Sikorsky experience indicate that the test 
data are in error by as much as 5°, so the stabilator control system was used to 
determine stabilator position for all static-trim comparisons. This method was also 
used in the Sikorsky validation (ref. 1). 

Level-flight trims were determined for airspeeds from 0 to 160 knots in 10-knot 
increments. The control positions, aircraft attitude, stabilator position, and 
main-rotor required power are presented in figure B1. The trim points from the 
Sikorsky simulation were provided by Sikorsky for comparison. Agreement between the 
two simulations is excellent for all airspeeds. As indicated in reference 1, agree- 
ment between the simulation and the flight test is good except for a nearly constant 
pedal bias of approximately 0.5 in. to the right. There is less agreement at lower 
speeds, caused in part by the pilot's difficulty in achieving trimmed flight on the 
back side of the power-required curve and also by the simplified modeling of the 
impingement of the main rotor downwash on the fuselage and stabilator. This also 
results in an underestimation of collective stick trim position by both simulations 
when in hover. 

The climb and descent comparisons are shown in figures B2 through B4 for air- 
speeds of 60, 100, and 137.5 knots. The Sikorsky simulation data are from refer- 
ence 1. Discrepancies exist between Ames and Sikorsky simulations because Sikorsky 
used a level-flight trim value of stabilator position for each speed. Again, agree- 
ment is good, with the greatest differences occurring in pedal positions. Figure B5 
contains comparisons for steady high-g coordinated turns for 60 and 100 knots. 
Because Sikorsky positioned the stabilator based on level-flight trim, differences 
between the simulations exist in pitch attitude and longitudinal cyclic position. 


DYNAMIC RESPONSE 


For transient response comparisons, the model configuration and test conditions 
were set up to duplicate those of the flight tests. All stability augmentation and 
flightpath stabilization was disabled and the pitch bias actuator was disabled and 
centered. As was done for static trim comparisons, the stabilator control system 
was used during the program trimming sequence to determine stabilator position, 
which was then held fixed. Because no information concerning moments and products 
of inertia of the test aircraft was available, the inertia terms used were those of 
the Sikorsky validation; they represent those of a typically configured and loaded 
Black Hawk. The ambient pressure, density, and temperature were determined from 
density and pressure-altitude test data. 

After trimming to the test conditions, time histories of test-aircraft control 
inputs were used as direct input to the simulation. Constant biases were added to 
the control inputs to correct for small discrepancies between the test aircraft and 
simulation trim positions. In this way, the simulation begins each transient 
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response in trim and simulation trim errors have a minimal influence on dynamic 
response comparisons. The test data as received were already processed and no 
additional adjustments were performed. In a few of the test time-histories, 
spurious data and obviously incorrect offsets can be seen. They are attributed to 
instrument biases or recording errors, but no effort has been made at correction. 

All simulation runs performed by Sikorsky in their nonreal-time validation used 
a time step of 0.01 sec. The simulation output was recorded only in time steps of 
0.03 to 0.05 sec because of plotting limitations as described in reference 1. The 
M/rev frequency characteristic of the simulation is higher than the Myquist cri- 
terion frequency corresponding to these sample rates, resulting in aliased output. 
The aliased frequencies seen in the Sikorsky simulation plots are 2.8 and 7.8 Hz and 
are especially prevalent in acceleration plots. Because of this, Sikorsky chose to 
base their simulation fidelity assessments on rate- and attitude-response data 
(ref. 1). A complete documentation of flight-test transient responses and Sikorsky 
simulation comparisons may be found in reference 11. 

Appendix C contains aircraft configuration and sensor location information. 
Appendix D presents transient-response time-history comparisons between flight test 
data and the Ames real-time simulation. The flight-test velocity and acceleration 
time-histories represent values measured at a sensor location on the aircraft. 
Simulation velocities and accelerations were computed at these locations so that a 
direct comparison can be made. It should be noted that blade lag- and flap-angle 
time histories represent the motion of one instrumented rotor blade. Because there 
is no direct correspondence in blade azimuth between simulation and flight test, 
only relative amplitudes and frequencies should be compared. 


Response to Control Inputs 

The transient responses contained in appendix D serve as examples of response 
to control inputs for hover, 60-knot, and 100-knot initial level-flight condi- 
tions. Initial responses compare well with flight-test data, with somewhat less 
damping of simulation rates in general. Aerodynamic damping of the fuselage is not 
modeled in the simulation j this may have some effect on vehicle response involving 
large angular rates. Because of pilot difficulty in maintaining trim of the 
unaugmented aircraft, in many cases flight-test data had drifted from trim at the 
time of the control input, causing differences between test and simulation 
responses. There is also evidence of control input hysteresis in flight test 
data. Small inputs sometimes produced no flight test response. Because there is no 
modeling of control freeplay, simulation response may be expected to be faster with 
greater magnitude. The simulation also produces off-axis responses of greater 
magnitude than flight-test responses, caused at least partly by model approximations 
to the mechanical control linkages. Control mixing is approximated as linear cross- 
couplings between control inputs in the simulation, whereas in fact the mechanical 
linkages are not linear. A list of model deficiencies may be found in reference 1, 
most of which are considered to be of minimal importance. Some of the deficiencies 


listed have been eliminated by model expansion and improvements, as described in the 
following section. 


Effects of Model Improvements 

Appendix E contains example time-history comparisons between the Ames real-time 
simulation and the Sikorsky simulation. Figure El, showing the response to a 
0.5-in. aft cyclic step in hover, illustrates the very good agreement between the 
two simulations under circumstances in which model updates have a small influence. 
The trim control positions are identical for both simulations and the velocity, 
position, and attitude responses are very similar. The small differences between 
the two simulations are caused by the modification to the primary servo dynamics. 

As explained previously, aliasing of output is evident in Sikorsky acceleration 
plots and amplitudes of oscillation differ between the two simulations because of 
differing coding implementations. The control input caused little variation of 
rotor speed, so blade motions and directional axis responses are in agreement. 

The response to a 1-in. down collective input at hover illustrates the improve- 
ments caused by program expansion and modifications (fig. E2). Modification to the 
primary servo bandwidth causes large improvements in all acceleration plots. As is 
most clearly seen in vertical-acceleration time histories, the Sikorsky simulation 
response lags the flight test response by approximately 0.1 sec, resulting in a 
simulation which causes pilot induced oscillation during hover. Increasing the 
servo bandwidth results in a very good agreement between the Ames simulation and 
flight-test data. 

Several of the simulation deficiencies described in reference 1 are a direct 
result of insufficient detail in engine- and drive-train modeling. Improvements to 
the Ames real-time model are shown by figure E2 in rotor speed, engine-related time 
histories, and yaw-related histories. The Sikorsky-simulation rotor speed varies 
adversely in response to a collective input because of an insufficient level of 
sophistication of the partial derivative engine model and fuel control system 
(detailed in ref. 2). Its response is also too fast and is of the wrong magni- 
tude. These rotor speed errors manifest themselves in yaw-acceleration, rate-, and 
heading-time histories as well as altering blade-lag response. Both engine models 
appear to trim at an incorrect level of power when compared with the flight test. 

This is a direct result of error in modeling of the main rotor torque in hover, as 
previously described. However, the NASA-developed engine shows slightly better 
agreement because of the more detailed representation of power required in the Ames 
real-time model. The simplified Sikorsky model always estimates tail-rotor required 
power to be 10% of main rotor power. In this hover case, tail rotor power of the 
real-time simulation is determined to be 13% of the main rotor power, which is a 
better approximation of the actual aircraft. 

An example of response at high speed is shown in figure E3. Again, accelera- 
tion and rate responses to the lateral input are lagged in the Sikorsky simulation 
whereas those of the Ames simulation are in good agreement with test data. Ames 
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simulation engine trim is in slightly better agreement with test results, as evi- 
denced by gas generator speed. The figure also illustrates the effect of slightly 
differing stabilator incidence angles on aircraft trim. The Sikorsky validation 
used the recorded flight test value of stabilator position for all transient 
response comparisons. Because the recorded stabilator data are suspected to be 
incorrect, the Ames validation was performed with an incidence angle determined by 
the stabilator control system during program trim. Blade-flap- and lag-time- 
histories, as well as trim control positions, suggest that this is a better approxi- 
mation to stabilator position. In this example, the control system positioned the 
stabilator at 3.75° trailing edge down, while the flight test data indicated a value 
of 7.8°. This caused a larger nose-down trim attitude in the Sikorsky simulation, 
resulting in lower magnitudes of first-harmonic flapping and lagging of the main 
rotor. The conclusion in reference 1 that lagging amplitudes are too small is 
actually the result of differing stabilator positions between simulation and flight- 
test data. 


RECOMMENDATIONS FOR IMPROVEMENT 


Several deficiencies noted in reference 1 have been eliminated by the more 
detailed representation of the rotor speed degree-of-freedom and modification to the 
linear representation of the primary servo dynamics. Those eliminated include the 
incorrect fuel flow, main rotor torque, and main-rotor-speed response to a collec- 
tive input, incorrect blade-flap and lag response to rotor control inputs, adverse 
yaw response to collective inputs, and delayed vehicle response to all control 
inputs . 

The minor off-axis couplings may be caused in part by a linear representation 
of the mechanical control linkages. The couplings are reduced by the stability 
augmentation system; they produced very few pilot complaints during a recent VMS 
simulation. Differences in trim pedal position and trim collective position in 
hover probably have little adverse effect, but they do point to areas where 
increased model detail may be desired. 

However, validation results and pilot opinion suggest that there is some need 
for further refinement and expansion of the model. Reference 1 states that large 
control inputs produce exaggerated model responses. Control freeplay should be 
modeled and a more sophisticated nonlinear representation of primary servo response, 
including rate and travel limits for each servo, may be warranted. This requires 
addition of the geometry that relates primary servo positions to swashplate 
angles. Pilot comments also suggest the need for simulation of the increased N/rev 
vibration level during quick-stops and in the speed range of 18 to 25 knots. This 
has been reported as a shortcoming of the actual aircraft (ref. 10), but pilots have 
found it to be an important cue, especially in the reduced visual environment of 
flight simulators. The simplified model of main rotor dynamic inflow may prevent 
physical modeling of the phenomenon, but an adequate representation may be achieved 
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through a flight-test-derived modification factor, as was done to simulate ground 
effect. 


CONCLUSIONS 


Real-time simulation of rotorcraft is achievable at the level of sophistication 
of the UH-60A version of Gen Hel, although a high-speed computational capability is 
required. Blade-element rotor models with rigid blades which are used in conjunc- 
tion with rigid-body airframe models may be executed in real time with no loss of 
model fidelity. Considerable effort is required to solve procedural and numerical 
problems which arise from the necessity to execute in the most computationally 
efficient manner possible. Because of the increasing use of high-fidelity modular 
programs, some form of generalized methodology is needed to decrease the development 
effort for each new program. Documentation of special correction techniques, such 
as that applied to the blade-lag-acceleration term, is necessary. Because it is 
generally not possible to achieve temporally perfect one-pass coding, the areas 
which must have correct temporal indices should be documented and incorporated into 
a generalized real-time modular program structure. 
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APPENDIX A 


MODEL UPDATES 


Modifications to the manufacturer-supplied mathematical model (ref. 2) are 
presented below. Symbols and program variable mnemonics follow the convention of 
the reference. 


ADDITIONS TO THE FLAPPING ACCELERATION EQUATION 

Additions to reference 2 are underlined. 

- M b , 

B = T“ <cos 3 [ V + e{2a(p cos ip - q sin iji) + p sin tp + q cos ip}] 

\ s s s s 
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^ s s s 

- sin <$[2(a + £)(p sin ip + q cos >p) + q sin ip - p cos ip ] } 

S S S S 


- sin e cos 8[(r - a ) 2 - 26(r - a) + <$ 2 ] + (M 


FAB 


M + M ) 
FFD FLD ; 


where 

8 blade flapping angle, rad 

6 blade lagging angle, rad 

e blade hinge offset from center of rotation, ft 
1^ inertia of blade about the hinge, slug*ft^ 

M b first mass moment of blade about the hinge, slug -ft 


( A 1 ) 
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^FAB aer °dynainic moment about hinge (blade span axis), ft ‘lb 
MppD flapping moment caused by flap damper, ft -lb 
MpLD flapping moment caused by lag damper, ft -lb 
8 rotor shaft speed, rad/sec 

^ rotor blade azimuth position, deg 

p s angular rate about shaft x-axis, rad/sec 

q s angular rate about shaft y-axis, rad/sec 

r angular rate about shaft z-axis, rad/sec 

V v shaft x-axis velocity, ft/sec 

x s 

V shaft y-axis velocity, ft/sec 
J's 

V shaft z-axis velocity, ft/sec 
z s 

Engine and Fuel Control Models 

Real-time General Electric T700 engine and fuel control models are described in 
Ballin, Mark G., "A High Fidelity Real-Time Simulation of a Small Turboshaft Engine" 
(in preparation). 


Drive Train Model 

The following description of the drive train and tail rotor torque computation 
is modified from unpublished Sikorsky documentation (D. G. C. Ruttledge, personal 
communication, Sikorsky Aircraft, United Technologies Corp., 1986). Simplifications 

made for real-time implementation include rigid drive shafts and a single-engine and 

clutch representation of the propulsion system. 

The drive train model consists of an idealized gearbox which sums the torque 
loads of a variety of components. Torque is provided by the engine module through a 

clutch interface. The gearbox model applies its torque load to the body c.g. via 

the rotor force and moment transformations rather than directly. All component 
loads other than internal torque balances are handled in the same manner as before. 

Gearbox - The idealized gearbox is represented by a large ring-gear which has 
inertia and oil-drag damping. All component torques are summed at this ring-gear. 

A body-centered axis set is employed — thus body yaw acceleration (r) terms are seen 
in the equations. 
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The ring gear inertia represents all the gear trains actually employed in the 
system, referred to a speed OMGBX which is the main-rotor-output shaft speed 
(rad/sec). To preserve the modularity concept of Gen Hel, components are specified 
with inertias including their associated shafting up to (an arbitrary) junction with 
the gearbox. This inertia is referred to the component operational speed. A gear 
ratio (n) between the component operating and the gearbox main-rotor output speed is 
also specified. The gearbox algorithm then sums the n 2 J appropriately, to repre- 
sent the component inertia. If a component contains a clutch mechanism (e.g., 
engines usually have overrunning clutches) the module representing that component 
must have the ability to add the additional inertia when the clutch is fully engaged 
and conversely to remove it as soon as the clutch begins to slip. 

In initial-condition mode: 

If OMGIBX = 0, OMGBX = OMRMR-OMGTMR 

Else OMGBX = OMGIBX 

In any case OMRBX = 0MGBX/0MGTMR 

In operate mode: 

no . of 
engines 

SQENBX = QQENG 

i = 1 

no . of 
loads 

QGBX = ^ QQ k - SQENBX + OMGBX • DAMPBX 
k= 1 

no . of 
loads 

FACIBX = FAC k 

k- 1 


OMGDBX 


RS .MR - QGBX/JGBX 
1 + FACIBX/JGBX 


OMGBX = 0MGBX o + J OMGDBX dt 
OMRBX = OMGBX /OMGTMR 
HPDBX = DAMPBX -0MGBX 2 /550.0 


While trimming, use OMGDBX as a criterion. 


(A2) 

(A3) 
( A4) 
(A5) 
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Note. Torque load modules called by GEARBOX must contain: 


1. 

a). = n t • OMGDBX 

(A6) 

2. 

oj. = n. • OMGBX 

l l 

(A7) 

3. 

QQ. = n. • QLOAD. 

ill 

( A8) 

4. 

FAC. = J. • n 2 

ill 

(A9) 


Where 

DAMPBX viscous damping component of the gearbox caused by oil drag and other fac- 
tors, ft»lb*sec 

FAC^ inertia of the ith component with reference to the gearbox speed 
FACIBX partial sum of the gearbox inertias 
HPDBX horsepower absorbed in the gearbox, hp 

JGBX the rotating inertia of the gearbox components referred to the output speed 
OMGBX (defined later), slug* ft 2 

Jj rotating inertia of the jth component referred to the component operating 

speed, slug -ft 2 

n- gearing ratio between component-j operating speed and gearbox output speed 

OMGBX 

OMGBX speed of the ring gear (output speed of main rotor gearbox), rad/sec 
OMGDBX acceleration of the ring-gear, rad/sec 2 
OMGTMR trimmed main-rotor speed (defined as input), rad/sec 
OMRMR main-rotor-speed ratio (defined with respect to OMGTMR) 

OMRBX gearbox-output-speed ratio ( defined . with respect to OMGTMR) 

QLOADj^ torque load of the ith component at the component operating speed, ft -lb 

QQ^ torque load of the ith component on the ring gear, ft*lb 

QQENG output torque of an individual engine (measured at the ring-gear), ft-lb. 

QGBX total torque load on the ring gear, ft* lb 
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p 

RS.MR body yaw acceleration, rad/sec 

SQENBX sura of the engine torques transmitted to the gearbox, ft ♦lb 

Tail Rotor Torque Transfer 

The tail rotor is treated as a typical torque-load submodule except that the 
gear ratio is not entered directly; rather it is computed from the nominal-trim 
rotor rotational speeds. 

In initial condition mode: 

n TR = OMEGTR/OMGTMR (A 10) 

In operate mode: 

QQTR = n TR * QTR 
QGBX - QGBX + QQTR 
FAC tr = n^ R * JSHTR 

Note that all forces and moments (including the pitching moment due to QTR) are 
transferred to the c.g. by the tail-rotor module itself. 

Where 

FAC tr inertia of the tail rotor with reference to the gearbox speed, slug -ft 2 

JSHTR inertia of the tail rotor and associated shafts and gears as reflected at 

the tail rotor hub, slug -ft 2 

n TR - overall gear ratio between tail rotor and gearbox main rotor output 
OMEGTR nominal tail-rotor rotational speed (defined as input), rad/sec 
OMGTMR nominal main-rotor trimmed speed (defined as input), rad/sec 
QGBX torque load of the gearbox ring-gear, ft -lb 

QQTR torque load of the tail rotor, reflected at the ring gear, ft -lb 

QTR torque load of the tail rotor, at the tail rotor hub, ft -lb 
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Main Rotor Torque Transfer 


The main rotor is treated as a typical torque-load submodule except that the 
gearing is assumed to be unity. 

OMR. MR = OMRDBX 

OMRMR = OMRBX 

FAC mr = JSHMR 

QQMR = QHMR 

QGBX = QGBX + QQMR 

Additionally, QHEG = QQMR (this allows the transfer of QQMR to the body c.g. in 
rotor code) . 

Where 

FAC mr inertia of the main rotor shafting and hub with reference to the gearbox 
speed, slug -ft 2 

JSHMR inertia of the main rotor shafting and hub, slug -ft 2 
OMR. MR nondimensional acceleration of the main rotor 

OMRDBX nondimensional acceleration of the main rotor, gearbox output, 1/sec 2 

OMRMR main-rotor-speed ratio (defined with respect to OMGTMR ) 

OMRBX gearbox-output-speed ratio (defined with respect to OMGTMR) 

QQMR torque load of the main rotor at the ring-gear, ft -lb 

QHMR torque load of the main rotor at the hub, ft -lb 

QGBX total torque-load at the gearbox ring gear, ft -lb 

QHEG yawing moment applied to the c.g. caused by the main rotor, ft- lb 

Accessory loads - The accessory load is a rather arbitrary torque load which 
represents an average load caused by accessories such as generators and hydraulic 
pumps. The nominal accessory hp is converted to a torque at trimmed rotor speed and 
a variation of 10% load per 100% rotor-speed change is assumed. 
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In initial condition mode: 


QQTACC = ACCHP * 550.0/OMGBX (All) 

In operate mode: 

QQACC = QQTACC (1.0 + 0.1 (OMRBX - 1.0)) (A12) 

QGBX = QGBX + QQACC 

FflC flCC = J ACC * n ACC 

where 


ACCHP 
FAC 


ACC 


ACC 


n 


ACC 


OAMGBX 

OMRBX 

QGBX 

QQACC 

QQTACC 


nominal accessory load in trim, hp 

inertia of the accessories with reference to the gearbox speed, slug*ft 2 
inertia of the accessories referred to the accessory drive speed, slug -ft 2 
gear ratio between accessory drive speed and gearbox output speed 
gearbox output speed, rad/sec 
nondimensional gearbox-output speed 

cumulative total-torque load at the gearbox ring gear, ft* lb 

actual accessory torque-load at gearbox ring gear, ft -lb 

torque load at the gearbox ring gear caused by accessories in nominal trim 
conditions, ft -lb 


Multiple engine clutch - This is modeled as a triple-channel overrunning clutch 
which allows the hub to freewheel at higher speed than the engine-shaft-output speed 
while never allowing gearbox speed to be greater than hub speed. 

The speeds are checked to see if they are in close proximity. If they are 
close, the accelerations are checked to see that the engine is applying a driving 
torque to the gearbox. Clutch chatter is reduced by delaying QPTE by one interval 
so that there are no mixed temporal indices in the clutch-logic equations. 

The clutch logic is: 


If OSEG < Mg * OMGBX, the clutch is disengaged 

Else, if Mg * OMDBX > QPTE/JE, the clutch is disengaged 

Else, the clutch is engaged 
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Whereupon 


If the clutch is engaged, 


FAC e = n* J E 
OSEG = OMGBX-n, 


OS. EG = OMGDBX *ng 
QQSTO = QPTE-n E 

but if the clutch is disengaged, FACIBX is not changed and 

OS. EG = QPTE/JE 


OSEG = OSEG, 


QQSTO r 0 


; 


OS. EG dt 


( A 1 4 ) 
( A 1 5 ) 
(A16) 
( A 1 7 ) 

( A 1 8 ) 
( A 1 9 ) 


In either case, QGBX = QGBX + QQSTO. 

All of the steps are repeated for the number of engines required or the QGBX 
equation is modified to: 

QGBX = QGBX + ENCOMP * QQSTO 
Also, if the clutch is engaged, 

FAC e = (ENCOMP - 1.0) • J E (A20) 

and the clutch logic is executed only once. 

ENCOMP number of identical engine/clutch combinations in simulation 

FAC E inertia of engine and associated shafting with reference to the gearbox 
speed, slug -ft 2 

J E engine and associated shafting inertia, slug*ft^ 

n E gear ratio between engine output and gearbox ring gear 

OMGBX gearbox-ring-gear-output speed 
OMGDBX acceleration of the ring gear 
OSEG engine-shaft-output speed, rad/sec 
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OS. EG engine-shaft-output acceleration, rad/sec^ 

QGBX cumulative total of ring-gear torques, ft -lb 

QPTE engine-output torque (referred to engine speed), ft -lb 

QQSTO engine torque (referred to the ring gear speed), ft-lb 


Drive Train Constants 

The drive train constants that were used for the UH-60A Black Hawk are as 
follows. 

DAMPBX =50.0 
JGBX = 208.4 
JSHMR = 41.85 
JSHTR = 19.69 
OMEGTR = 124.62 
ACCHP =40.0 

J ACC = °- 00 ^ 

n ACC = 51 
n E = 81.042 

JE = 0.124 
ENCOMP = 2 


ADDITIONS TO TAIL ROTOR MODEL 
(Modified from unpublished Sikorsky documentation) 


The tail-rotor torque coefficient is defined as 
CQTR = - t ( a 2 t 5 5 - a TR t 41 )X + {6 2 t 5.6 " a TR t 4.2 )8 TR X + (S 2 fc 5.7 ' a TR t 4.3 )6 1 X 

* (5 2 t 5.8 - a TR t 4.4 )0 TR 2 + (S 2 fc 5.9 * a TR fc 4 ,5 ] 0 1 9 TR 


+ (s 2 fc 5.10 ‘ a TR fc 4.6 )0 1 + S oS.1 + 6 1 t 5.2 X + 6 1 t 5.3 9 TR + 6 1 fc 5 ,4 0 1 
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where a = bc/irR^ and 0-j = TWSTTR/57 .. 3 . Tail rotor torque is then computed as 

QTR = CQTR-upn^R^OMGBX 2 

The coefficients used in the above equations are defined as 

2 2 2 
ET 5 y B* 2 

C 4. 1 = 2 + 4 + 1296 U TR 


t - — 

C 4.2 " 3 


2 Q 

8B Y B* 2 
3 + 864 U TR 


'4.3 


= r- + 2B 


2 y 2 B 1Q 2 
+ 1080 V TR 


t 8 „2 y 2 B 1Q 2 

C 4.4 '9 + 2304 U TR 


4B 2 yV 1 1 2 


'4.5 = 3 + 1440 m TR 


'4.6 


4 2 12 

B_ y B 2 

2 + 3600 U TR 


-B 5 2 

'4.7 = 108 U TR 


-B 6 2 

'4.8 " 144 y TR 


7 

-B' 2 

'4.9 = 180 m TR 


'4.10 


2 

B_ 2 
36 U TR 


1 1 2 

'5 . 1 = 4 + 4 m TR 
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fc — — 

5.2~ 3 

. 112 

5.3 * 4 + 4 y TR 


t 


5.4 


1 1 2 

5 + 6 y TR 


4 5 6 

. i zl L_ J 2 B B B 2 

5.5 = 2 + 4 + B 2 + 2B 4 + Y 162 " 81 + 144 U TR 

2 4_ 4 2 B^_ B^ B^ 2 

t 5.6 = 3 + 3B * D 3 + Y 108 “ 54 + 96 U TR 

*. 1-1 2 B 6 2B 7 B 8 2 

0.7 ~ 2 + 1 + b 2 + Y 135 ” 135 + 120 y TR 

1 i q p r 6 r 7 r 8 ~ 

fc 5.8 = 4 + 4 + T72 +y 288 “ T44 + 256 U TR 

yt5 

. 2 1 1_ 2 B 7 B 8 B 9 2 

5.9 " 5 + 3 + 3B + Y 180 - 90 + 160 M TR 

15 2 B 8 B 9 B 10 2 

5.10 = 6 + 8 + Y 450 " 225 + 400 U TR 


for which 

Y = 

where 

a^ R blade-section lift curve slope, 1/rad 

b number of tail rotor blades 

B blade-tip-loss factor 

c tail-rotor-blade chord, ft 


pa TR C TR R TR 


IB 


TR 
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CDTR 


p 

tail rotor hub drag, ft 

p 

IB tr mass moment of inertia of rotor blade, slug* ft 

OMGBX gearbox ring gear output speed, rad/sec 

QTR tail rotor torque, ft*lb 

R-p^ tail rotor radius, ft 

TWSTTR linear blade twist, deg 

5 0 ,l5 1 ,l5 2 coefficients in power series which express section profile-drag coefficient 
as a function of lift coefficient 

Y mass constant of tail rotor blade 

X tail rotor inflow 

tail-rotor-tip-speed ratio 

S2-PP tail-rotor trim speed, rad/sec 

p air density, slugs/ft J 

a tail rotor solidity 

0 1 linear blade twist, rad 

ADDITIONAL TAIL ROTOR CONSTANTS FOR UH-60A BLACK HAWK 


The constants necessary to incorporate the tail rotor torque model are as 
follows : 


5 = 0.0087 

o 

6 = - 0.0216 
S 2 = 0.4 

IB tr = 3.1 


In addition, the value of CDTR in reference 2 has been changed: 


CDTR = 0.4 
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PRIMARY SERVO DYNAMICS MODIFICATION 


The transfer function as it appears on pp. 5.5-14, 5.5-15, 5.5-16, and 5.5-17 
of reference 2 is as follows. 

S MX 1 

6 MIX 0.00114 s 2 + 0.0473 s + 1 

The primary servo dynamics modification as based on no-load acceptance testing 
of actuators: 

5 MX _ 1 

6 MIX 0.0002041 s 2 + 0.020 s + 1 
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APPENDIX B 


STATIC TRIM CHARACTERISTICS 


Steady-stafce trim characteristics of the Ames and Sikorsky versions of Gen Hei 
are compared with USAAEFA test data contained in reference 10. 


GROSS WEIGHT 160001b 

CG STATION 351.0 in. 

CG WATERLINE 231.5 in. 

DENSITY ALTITUDE 5250 ft 


SIKORSKY GEN HEL 

AMES GEN HEL 

O AEFA TEST 



100 

80 


r rt 



_ RT 

o Q- 

o 

o 


(e) 

1 




o 

i 


o 



for 

,NGLE, deg 


TRAILING EDGE DOWN 






% *>IDI1S 3AI1D3TTOD 


% 'SlVQ3d 


















60 knots 



100 knots 

FLT16 RUNS 41-45 






Figure B5.- UH-60A Black Hawk steady coordinated-turn trim, 60 knots and 100 knots. 
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Figure B5.- Concluded. 
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APPENDIX C 


TEST CONDITIONS FOR TRANSIENT RESPONSES 


Test conditions and aircraft configuration information are presented below for 
each of the transient-response time histories of appendices D and E. Pertinent 
inputs to the simulations are also given when they differ from or are estimates of 
flight-test values. 


SENSOR LOCATIONS 


CG Accelerometer 

Station. .389 in. 

Buttline .-31 in . 

Water 1 ine 207 .7 in . 

Translational Velocity Transducer 

Station. .248 in. 

Buttline.. 70 in. 

Waterline 265 in. 


TEST DATA CONDITIONING 

Samples per second 100 

Filter break frequency. 30 Hz 


TEST CONFIGURATION FOR ALL TRANSIENT RESPONSE COMPUTATIONS 


, I xx = 4659 
Iyy = 38512 
I zz = 36796 
*xz = 1882 


slug -ft 2 
slug -ft 2 
slug* ft 2 
slug -ft 2 


Buttline CG location 
Pitch bias actuator 
Flightpath stabilization 
Stability augmentation system 
Stabilator position 


0 in. 

disabled and centered 

disabled 

disabled 

fixed to the trim value set by control system 
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TABLE Cl.- AIRCRAFT CONFIGURATIONS AND TEST CONDITIONS 
FOR TRANSIENT RESPONSES 


File Flight 

Run 

Gross 

Weight, 

lb 

CG 

Station, 

in. 

CG 

Waterline, 

in 

Density 

Altitude, 

ft 

Pressure 

Altitude, 

ft 

Equivalent 

Airspeed, 

knots 

Rotor 

Speed, 

rad/sec 

102 

18 

20 

15990 

348.8 

231.7 

6809 

5500 

93.0 

27.0 

103 

18 

10 

16150 

349.5 

231.6 

6565 

5256 

91.9 

26.86 

201 

47 

27 

15870 

359.3 

232.6 

6658 

5582 

0 

27.0 

202 

47 

31 

15760 

359.2 

232.7 

6673 

5493 

0 

27.0 

204 

48 

21 

15950 

359.4 

232.5 

6510 

5370 

0 

27.0 

205 

48 

23 

15870 

359.4 

232.6 

6777 

5584 

0 

27.0 

206 

48 

26 

15720 

359.2 

232.7 

6834 

5562 

0 

27.0 

209 

49 

35 

15750 

359.3 

232.7 

6659 

5742 

0 

26.93 

210 

50 

27 

15940 

359.4 

243.8 

6705 

5776 

0 

26.90 

212 

50 

24 

16030 

359.4 

232.5 

6296 

5411 

0 

26.87 

401 

45 

16 

16530 

359.6 

232.3 

5084 

4484 

93.3 

27.0 

403 

45 

44 

15840 

359.5 

232.6 

6270 

5641 

91.0 

27.05 

404 

46 

12 

16540 

359.5 

232.3 

5001 

4455 

92.8 

27.0 

405 

46 

14 

16490 

359.5 

232.3 

4872 

4401 

55.8 

27.01 

409 

43 

15 

16410 

359.6 

232.3 

4360 

4556 

93.3 

26.53 

412 

49 

16 

16470 

359.7 

232.3 

4995 

4554 

92.8 

26.95 

502 

69 

18 

16110 

351.1 

232.4 

6268 

4969 

56.9 

27.0 

503 

69 

20 

16050 

351.1 

232.4 

6313 

5020 

54.6 

27.0 

504 

14 

23 

15770 

348.3 

231.7 

5718 

4514 

54.4 

26.57 

505 

14 

24 

15760 

348.3 

231.7 

5796 

4636 

55.0 

26.59 

509 

25 

47 

15510 

352.0 

232.0 

6969 

5953 

53.1 

26.85 

510 

29 

30 

15720 

351.4 

231.8 

5930 

5515 

53.96 

26.90 

512 

29 

36 

15600 

351.4 

231.9 

5631 

5948 

53.1 

26.90 

706 

17 

48 

15410 

351.9 

232.1 

7427 

6450 

92.8 

26.90 

710 

18 

33 

15740 

347.6 

231.8 

7015 

5636 

89.8 

27.03 ■ 

711 

' 18 

37 

15690 

347.4 

231.8 

6988 

5615 

90.3 

27.01 
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APPENDIX D 


DYNAMIC RESPONSE CHARACTERISTICS 


Transient responses of the Ames real-time simulation and USAAEFA flight test 
data of reference 11 are compared for small step-inputs of each control. No stabil- 
ity augmentation or flightpath stabilization was used. One primary control was 
input for each case while other controls were held fixed as long as possible. If a 
step input was not possible for reasons of instability or safety, a pulse input was 
used. 
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PEDALS, in. COLLECTIVE STICK, in. 
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LONGITUDINAL 

CYCLIC, in. PITCH ACCELERATION, 1 /sec 2 PITCH RATE, deg/sec PITCH ATTITUDE, deg 


GEN HEL (AMES VERSION) 


- - PEER TEST DRTfl 



Figure D1.- Continued. 
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LATERAL CYCLIC, in. ROLL ACCELERATION, 1 /sec^ ROLL RATE, deg/sec ROLL ATTITUDE, deg 


GEN HEL (RUES VERSION) 


REFfi TEST DRTfl 




TIME (SEC) 




Figure D1.- Continued. 
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YAW ACCELERATION, 1/sec 2 YAW RATE, deg/sec HEADING, deg 


GEN HEL (AMES VERSION) 


- - REEF! TEST DATA 



Figure D1.- Continued. 
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Figure D1.- Continued. 
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BLADE FLAP, deg BLADE LAG, deg ANGLE OF SIDESLIP, deg ANGLE OF ATTACK, deg 


GEN HEL l RUES VERSION) FIEFA TEST DATA 



TIME I SEC) 



TIME (SEC) 



TIME (SEC) 



TIME (SEC) 

Figure D1.- Continued. 
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VERTICAL ACCELERATION 

ROTOR SHAFT HORSEPOWER C G SENSOR ft/sec 2 LATERAL ACCELERATION LONGITUDINAL ACCELERATION 

"7 " f C.G. SENSOR, ft/sec 2 C.G. SENSOR, ft/sec 2 


GEN HEL (RUES VERSION) 


REFR TEST ORTR 



TIME (SEC) 



TIME (SEC) 



TIME (SEC) 



Figure D1.- Continued. 
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27.0 



Figure D1.- Concluded 










PEDALS, in. COLLECTIVE STICK, in. 


GEN HEL l AMES VERSION) 


AEFA TEST DATA 





o.o o.s 



Figure D2.- Comparison of Ames real-time simulation and flight-test data. 
(File 202: 1-in. left cyclic input at hover) 
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LONGITUDINAL 

CYCLIC, in. PITCH ACCELERATION, 1/sec 2 PITCH RATE, deg/sec PITCH ATTITUDE, deg 


GEN HEL (AMES VERSION) 


AEEA TEST DATA 





TIME (SEC) 


S.O | — i— »— T —t— r- i — ■- ■ ■ i ■ ■ ■ ' i ' < ■ i- 1 ■ " »■■■ ' i i ■ ■ - i — i ' ■ ' ■ t ■ ' t— r-r ' ■ 1 ■ ■ ■ ' i 


2.5 L 


0.0 


-2.5 I- 


■5.0 1 • ■ ■ ■ 1 ■ ■ ■ .... i .... i ■ ... ■ .... i .... i . ... i .... i .... i .... I 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 <.5 5.0 5.5 6.0 6.S 

TIME (SEC) 


Figure D2.~ Continued. 
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LATERAL CYCLIC, in. ROLL ACCELERATION, 1/sec: ROLL RATE, deg/sec ROLL ATTITUDE, deg 


GEN HEL (RMES VERSION) 


REFR TEST ORTfl 



TIME I SEC) 



TIME (SEC) 




TIME (SEC) 


Figure D2.- Continued. 
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PEDALS, in. YAW ACCELERATION, 1/sec 2 YAW RATE, deg/sec HEADING, deg 


GEN HEL (RUES VERSION! 


AEFft TEST OflTfl 



TIME (SEC) 



TIME (SEC) 



TIME (SEC) 


.0 2.5 

3.0 

TIME 

3.5 4.0 4.S 

(SEC) 

s.o 

5.5 

6.0 

6.5 

Figure 

D2.- 

Continued. 
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Figure D2.- Continued. 


COLLECTIVE STICK, in 


VERTICAL VELOCITY 
SENSOR, ft/sec 



62.5 


LATERAL VELOCITY 
SENSOR, ft/sec 


LONGITUDINAL VELOCITY 
SENSOR, ft/sec 



GEN HEL (AMES VERSION) REFR TEST DRTR 






BLADE FLAP, deg BLADE LAG, deg ANGLE OF SIDESLIP, deg ANGLE OF ATTACK, deg 


SEN HEL l RUES VERSION) 


REFR TEST ORTR 



TIME (SEC) 



TIME (SEC) 



TIME (SEC) 



Figure D2.- Continued. 
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• GEN HEL l FINES VERSION! 



Figure D2.- Concluded. 
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GEN HEL (AMES VERSION) 
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Figure D3.- Comparison of Ames real-time simulation and flight-test data. 
(File 204: 0.5-in. forward cyclic input at hover) 
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Figure D7.- Comparison of Ames real-time simulation and flight-test data. 
(File 509: 1-in. aft cyclic input, 60 knots) 
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Figure D9.- Comparison of Ames real-time simulation and flight-test data. 
(File 510: 0.5-in. up collective input, 60 knots) 
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Figure DIO.- Comparison of Ames real-time simulation and flight-test data. 
(File 512: 1-in. down collective input, 60 knots) 
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Figure Dll.- Comparison of Ames real-time simulation and flight-test data. 
(File 502: 0.5-in. right pedal input, 60 knots) 
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Figure D12.- Comparison of Ames real-time simulation and flight-test data. 
(File 503: 1-in. left pedal input, 60 knots) 
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Figure D 1 3 . — 


Comparison of Ames real-time simulation and flight-test data. 
(File 401: 1-in. aft cyclic input, 100 knots) 
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Figure D 1 6 . - Comparison of Ames real-time simulation and flight-test data. 
(File 102: 1-in. right pedal input, 100 knots) 
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APPENDIX E 


SIKORSKY AND AMES GEN HEL DYNAMIC RESPONSE COMPARISONS 


Transient responses of the Ames real-time simulation, the Sikorsky Gen Hel 
simulation, and USAAEFA flight test data of reference 11 are compared for small 
step- inputs of each control. 
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Figure E2.- Comparison of Ames real-time simulation, Sikorsky simulation, and 
flight-test data. (File 210: 0.5-in. down collective input at hover) 
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Figure E2.- Continued. 
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Figure E2.- Continued. 
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Figure E2.- Continued. 
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Figure E2.- Continued. 
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Figure E2.- Concluded. 
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Figure E3.- Comparison of Ames real-time simulation, Sikorsky simulation, and 
flight-test data. (File 409: 1-in. right lateral input, 100 knots) 
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Figure E3-- Continued. 
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Figure E3.- Continued. 
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Figure E3.- Continued 
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Figure E3.- Continued. 
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Figure E3.- Concluded. 
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